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We consider the spin-1/2 two-leg ladders with ferromagnetic (FM) interactions along legs and 
rungs. Using the stochastic series expansion QMC method, we study the low-temperature magnetic 
behavior of the system. An isolated spin-1/2 FM two-leg ladder is in the gapped saturated FM phase 
at zero temperature. As soon as the spin-1/2 FM two-leg ladders are connected with antiferroniag- 
netic (AFM) inter-ladder interaction, a first-order commensurate-incommensurate quantum phase 
transition occurs in the ground state magnetic phase diagram. In fact a jump in the magnetization 
curve is observed. We found that, conpled spin-1/2 FM two-leg ladders are in a nonmagnetic phase 
at zero temperature. Applying a magnetic field, the ground state of conpled spin-1/2 FM two-leg 
ladders remains in the nonmagnetic phase up to a quantum saturate critical held. 

PACS numbers: 75.10.Jm; 75.10.Pq 


I. INTRODUCTION 

Recently, spin-1/2 two-leg ladder systems have de¬ 
voted considerable growth to itself experimentally and 
theoreticallyir— . The main interest to study of these 
systems is related to this fact that the high-Tc 
superconductivity^^— phenomenon occurs in these sys¬ 
tems and also they have a gap in the spin excitation 
spectrum. 

In the study of spin-1/2 two-leg ladder systems with 
antiferromagnetic (AFM) leg and rung interactions, the 
formation of spin singlets located on each rung open the 
spin gap in the energy spectrum which is called the gaped 
spin liquid phaseii^. These kind of AFM two-leg ladder 
systems are observed in the nature^r— . The effect of 
a magnetic field on the physical properties of these com¬ 
pounds has been a field of intense studies. It is found that 
there are two quantum critical fields in the ground state 
phase diagram of these kind of spin-1/2 two-leg ladders. 
Generally, at low magnetic fields (h < hc^), there is a 
spin liquid phase (a gapped phase) at low temperature^. 
Both magnetic susceptibility and the magnetization go 
up first with cooling, then decay exponentially to zero 
at low temperatures. Also, the specific heat has a single 
peak at low temperature due to transition from disor¬ 
dered phase to the spin singlet gapped phase— >2^. The 
Tomonaga-Luttinger liquid (TLL) gapless phase is found 
in hc„ < h < hc 2 regime at low temperatures^i^. One 
of the spin liquid {h < [hc^ hc^)/'^) or spin polarized 
{h > (hci + /ic 2 )/ 2 ) phases at higher temperature is ex¬ 
pected. The thermodynamic properties like magnetiza¬ 
tion and the susceptibility have a finite value at low tem¬ 
perature which show the vanishing of the energy gap in 
the TLL phase. Specific heat shows a second peak and 
goes down linearly with lowering temperature in the TLL 
regime^. 

The spin-1/2 two-leg ladder systems with AFM legs 
and FM rungs are also observed from experimental point 
of view22iiH. By means of the specific heat and the mag¬ 
netocaloric effect measurements, a phase boundary be¬ 


tween the spin liquid phase and the ordered phase is 
determined^!. We have to mention that, the ladders with 
FM leg and AFM rung exchange interactions found ex¬ 
perimentally during last two years22r— . 

Recently, Nagashiwa and coworker synthesized 
chlorido-bridged dinuclearcopper(II) complex with 
2-methylisothiazol-3(2H)-one with chemical formula 
Cl(0-mi)2(-Cl)2)}^- This compound is a quasi- 
two-leg ladder system with FM exchange interaction in 
legs and rungs. The Weiss temperature is estimated 
about 8.7 K, indicating ferromagnetic behavior. Also no 
exotic phenomena due to spin frustration were observed 
within the measured temperature range. From theoret¬ 
ical point of view, spin-1/2 two-leg ladders with FM 
legs and rungs are much less studied. In a very recent 
work, the temperature dependence of the magnetic 
susceptibility and the magnetic structure factors is 
studied using the modified spin wave theory and the 
numerical exact diagonalization technique!!. They have 
showed that in an intermediate temperature range, 
their analytical results are consistent with the numerical 
exact diagonalization results. The need for investigation 
of AFM inter-ladder coupling in FM ladders ( FM 
exchange interaction in both legs and rungs ) is also 
motivated by synthesizing the 3-C1-4-F-V, 3-Br-4-F-V 
and 3-TV crystals recently which are the candidates for 
two-leg ladders with FM legs and AFM rungs!!"— . 

In this paper we study the AFM inter-ladder interac¬ 
tion effect between FM two-leg ladders in the presence of 
a magnetic field (see Fig. I). This system for large inter¬ 
ladder interaction can be considered as the ladders with 
FM legs and AFM rungs by FM coupling. We used the 
recently developed stochastic series expansion (SSE)!!!! 
QMC method to provide numerical simulation results. 
The magnetization, the susceptibility and the specific 
heat are calculated for large enough finite size systems. 
Our simulation results show that isolated FM two-leg lad¬ 
ders are in an ordered phase at zero temperature. As soon 
as the AFM inter-ladder interaction is added, a first-order 
commensurate-incommensurate quantum phase transi- 
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FIG. 1: The schematic picture of FM two-leg ladders with 
AFM inter-ladder interaction. The size of coupled ladders is 
10 X 30. 

tion occurs from the gapped ferromagnetic phase to a 
nonmagnetic phase. By applying a uniform magnetic 
field, the ground state of coupled ladders remains in this 
nonmagnetic phase up to a saturated critical field. 

The outline of the paper is as follows: In section II we 
introduce our numerical simulation results for an isolated 
FM spin-1/2 two-leg ladder. In section III we consider 
coupled FM ladders and present results of our QMC sim¬ 
ulation on the low-temperature magnetic behavior of the 
system. Finally we conclude and summarize our results 
in section IV. 


II. AN ISOLATED FM TWO-LEG LADDER 

In this section we consider a FM spin-1/2 two-leg lad¬ 
der system in the presence of an external magnetic field. 
The Hamiltonian of the system is written as 

71,0. n,oc 

H” Jrung ^ ^ i 2 (1) 

n 

where „ is the spin 5 = 1/2 operator on rung n 
(n = 1,...,L) and leg a (a = 1,2). An applied mag¬ 
netic field h in the Z direction leads to zeeman term. 
The rung and leg exchanges are denoted by Jmng and 
Jieg, respectively. 

It is found^^ that the rung exchange Jmng = —9.68 K 
is about four times larger than Ji^g = —2.068 K. In the 
following, we consider these values in our QMC simula¬ 
tion approach. We have used the ALPS^ code which 
is known as one of the best codes in this field. We 
have considered two-leg ladders with the maximum size 



FIG. 2: The spontaneous magnetization as a function of the 
temperature for isolated spin-1/2 FM two-leg ladders. As it 
is seen, at zero temperature the ground state of the system 
is in a saturated FM phase. By increasing the temperature 
the induced thermal fluctuations decrease the magnetization 
of the system. In the inset, the magnetic susceptibility in 
the absence of the magnetic field is plotted as a function of 
the temperature. QMC simulation has been carried out for 
Heisenberg model. The size of two-leg ladder is 2 x 100. 


N = 2L = 2 X100 and periodic boundary conditions. The 
QMC simulation is performed for the maximum 1000000 
equilibration sweeps and 2000000 measurement steps. 

Fig- m presents the magnetization as a function of the 
temperature in the absence of the magnetic field. In¬ 
duced thermal fluctuations by increasing temperature 
cause in decreasing the magnetization and it reaches zero 
for high enough temperatures. The magnetic suscepti¬ 
bility is also shown in the inset of Fig. [2l It only goes 
up with cooling the system which is in good agreement 
with the experimental results^. The asymptotic behav¬ 
ior confirms that the ground state of the system is in 
a magnetic ordered phase at T = 0, which is also in 
complete agreement with our results on the magnetiza¬ 
tion. In addition, we have also studied the effect of the 
magnetic field on the thermodynamic behavior of an iso¬ 
lated FM two-leg ladder system. Our QMC simulation 
results are presented in Fig. [3Da) for different values of 
the magnetic field h = 0.5 T, 1.0 T, 2.0 T, and 3.0 T. 
we have calculated the spin gap according to expres¬ 
sion x(T) ^exp(—A/T)/-\/T'. The obtained results are 
A = 0.88 K, 1.6 K, 2.6 K, and 4.5 K, for h = 0.5 T, 
1.0 T, 2.0 T, and 3.0 T, respectively. 

As shown in Fig. [SDa), the magnetic susceptibility goes 
up first with lowering temperature until it reaches to a 
maximum value, then decreases exponentially down to 
zero at low temperatures. This exponential fall of sus¬ 
ceptibility indicates that there is a gap in the excitation 
spectrum of the FM two-leg ladder system in the pres¬ 
ence of a magnetic field. With decreasing the magnetic 
field the peak of susceptibility curve increases and thus 
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FIG. 3: (a) Magnetic susceptibility x[T) versus temperature 
for isolated spin-1/2 FM two-leg ladders at different values 
of the magnetic field. QMC simulation has been carried out 
for Heisenberg model, (b) Magnetic specific heat per site ver¬ 
sus the temperature for different values of the magnetic field. 
QMC simulation has been carried out for Heisenberg model. 
The size of two-leg ladder is 2 x 100. 


the spin gap is weakened. 

To complete our study on the thermodynamics of an iso¬ 
lated FM spin-1/2 two-leg ladder, we have also calculated 
the specific heat, Cm{T). Our numerical results are pre¬ 
sented in Fig. EJb) for different values of the magnetic 
field. As it is seen, at zero temperature the specific heat 
is zero, by increasing the temperature it remains almost 
zero up to a threshold temperature which is known as the 
indication of the spin gap. We have also checked the low 
temperature results and found an exponential behavior 
which is also an indication of the gapped phase. 


High-held nuclear magnetic resonance and inelastic 
neutron scattering measurements in some quasi-two-leg 
ladder system a^^d^i^5,26 weak coupling be¬ 

tween ladders induces ZD phase transition at low temper¬ 
ature region. In these compounds, phase transition be¬ 
tween quantum phases and the 3D gapless-Neel ordered 
phase occurs at very low temperatures. 

To the best of our knowledge, the effect of the inter¬ 
ladder coupling in spin-1/2 two-leg ladders with FM legs 
and rungs has not been studied so far. Interesting re¬ 
sults are not expected for the inter-ladder FM coupling. 
For this reason, we have considered the inter-ladder AFM 
coupling. We have performed the QMC simulations on 
the AFM coupled FM two-leg ladder systems for differ¬ 
ent values of the inter-ladder exchange Jin and in the 
presence of a magnetic held. The scheme of coupled lad¬ 
ders is illustrated in Fig. [TJ The numerical results are 
calculated for coupled ladders with Jmng = —9.68 AT, 
Jieg = -2.068 K, Jin/Jieg = -0.1,-0.2,...,-0.5 and 
different sizes up to 5 x 2L = 5 x (2 x 30). 

Fig. HD a) shows the magnetization, M, of the coupled 
system as a function of the magnetic held in the vicin¬ 
ity of the zero temperature and Jin/Jieg = —0.2. It is 
clearly seen that the magnetization is zero in the absence 
of the magnetic held which shows that the induced quan¬ 
tum huctuations by inter-ladder interaction destroy the 
long-range FM order of the ground state of isolated FM 
two-leg ladders. It means that in the ground state phase 
diagram of the coupled FM ladders, the line Jin/Jieg = 0 
must be a quantum critical line. By applying the mag¬ 
netic held, the magnetization starts to increase and will 
be saturated in a critical magnetic held, which depends 
on the value of the inter-ladder exchange interaction. 
So, such behavior of magnetization upon contribution of 
AFM inter-ladder interaction suggests, the existence of a 
nonmagnetic phase. 

Now let us see what happens when we attempt to 
change the inter-ladder interaction at zero temperature 
and in the absence of the magnetic held. Fig. HDb) 
presents the magnetization versus inter-ladder Jin in¬ 
teraction in the absence of the magnetic held. The 
inter-ladder interaction effects embedded in the cou¬ 
pled ladders drastically change the value of magnetiza¬ 
tion. There is a sharp drop of magnetization in the low 
inter-ladder coupling about O.OU/eg. This behavior is 
known as the main indication of a hrst-order phase tran¬ 
sition. Therefore it suggests the existence of a hrst-order 
commensurate-incommensurate quantum phase transi¬ 
tion in the ground state phase diagram of the model. 
Moreover, to conhrm the type of the mentioned quan¬ 
tum phase transition in our model, we have also imple¬ 
mented the Lanczos algorithm to hnd the ground state 
energy of the system. A very important indication of 
the hrst order phase transition is the discontinuity in the 
hrst derivative of the ground state energy at the quantum 
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FIG. 4: (a) The magnetization of spin-1/2 FM two-leg lad¬ 
ders with inter-ladder AFM interaction as a function of the 
magnetic field in the vicinity of the zero temperature and 
Jin/Jieg = —0.2. QMC simulation has been carried out for 
Heisenberg model. The size of coupled ladders is 10 x 30. (b) 
The magnetization versus inter-ladder Jin interaction in the 
absence of the magnetic field. The size of coupled ladders is 
10 X 30. The effect of AFM inter-ladder interaction causes 
the system undergoes a first-order phase transition. QMC 
calculation has been carried out for Heisenberg model. In the 
above inset, specific heat per site was plotted as a function of 
the temperature in the absence of the magnetic field by using 
QMC simulation. In the bottom inset, using numerical lanc- 
soz algorithm, the first derivative of energy per site versus Jin 
was plotted. 


critical point. Using the numerical Lanczos method, we 
have calculated the ground state energy for system sizes 
2 X (2 X L) = 16, 20 and plotted the first derivative of the 
ground state energy per site as a function of the inter¬ 
ladder interaction in the bottom inset of Fig. 4(b). The 
results show a clear discontinuity in the first derivative 
of ground state energy which is in agreement with our 
QMC results. To find out the effect of inter-ladder AFM 
interaction, we have also presented the QMC results for 
the specific heat. Cm, in the top inset of Fig. IUb). In 


FIG. 5: (a) Magnetization versus temperature for spin-1/2 
FM two-leg ladders coupled with AFM inter-ladder interac¬ 
tion at different values of the magnetic field. QMC simulation 
has been carried out for Heisenberg model, (b) the magnetic 
susceptibility versus the temperature for different values of 
the magnetic field above and below quantum critical field. 
QMC simulation has been carried out for Heisenberg model. 
The size of coupled ladders is 10 x 30. 


this figure, we observe a second peak at very low temper¬ 
atures which can be considered as an indication of the 
mentioned nonmagnetic phase. 

It is interesting to discuss the thermodynamic prop¬ 
erties of coupled ladders in the gapless phase Jin = 
—0.2 Jieg for different magnetic fields. In general there is 
a similarity between isolated FM ladders and the AFM 
coupled FM ladders for high magnetic fields. But, the 
differences can be noticed within the range of 0.0 < h < 
0.5 T. Fig. inia) shows the magnetization versus the 
temperature for h = 0.2 T, h = 0.3 T, h = 0.45 T, 
and h = 0.5 T. The magnetization is not saturated for 
h < 0.45 T at zero temperature. To find the critical mag¬ 
netic field he, we have performed the calculation for other 
magnetic fields within the range of 0 T < h < 0.5 T. 
There is no saturation in magnetization curve up to 
































5 



FIG. 6: Magnetic specific heat per site versus the temper¬ 
ature for spin-1/2 FM two-leg ladders coupled with AFM 
inter-ladder interaction in the region of the magnetic field 
h < he — 0.42T. In the inset, the magnetic specific heat 
per site is plotted for values of the held more than the quan¬ 
tum critical held. QMC simulation has been carried out for 
Heisenberg model. The size of coupled ladders is 10 x 30. 

he = 0.42 T which is in accordance with Fig. IHa). In the 
inset of Fig. [5ja) we have plotted the magnetization as a 
function of the temperature to show that there is not any 
plateau in the curve of the magnetization. As shown in 
Fig. [Sjb) the magnetic susceptibility is plotted for differ¬ 
ent magnetic fields in the ratio of Jin/Jieg = —0.2. There 
are indications that the FM two-leg ladder spin systems 
under this weak AFM inter-ladder interaction undergo 
to a nonmagnetic phase in the low magnetic fields. x{T) 
goes to finite values at low temperatures within the range 
of h < he- Also this behavior suggests that the coupled 
ladders should be in gapless phase, whereas the exponen¬ 
tial fall-off of x(r) above critical field is a signature of 
existence of gap in this system. 

In Fig. [6] the curve of Cm{T) at h < he = 0.42 T 
and h > he (in the inset) is depicted.lt can be clearly 
seen that below the quantum critical field, a second peak 
appears at low temperatures, whereas as illustrated in 
the inset, above the quantum critical field he, the second 
peak disappears. These results indicate that the system 
undergoes a cross-over from gapped phase to a nonmag¬ 
netic phase. Also this behavior suggests that the coupled 
ladders should be in gapless phase. Generally, the weak 
AFM inter-ladder interaction has the substantial effect 
at low temperatures within the range of low magnetic 


fields h<he = 0.42 T. 


IV. CONCLUSION 

In summary, we have calculated the thermodynamic 
properties of FM two-leg ladder like Cl(0-mi)2(- 

01 ) 2 )]. We have performed stochastic series expansion 
QMC using ALPS code to investigate spin-1/2 isolated 
FM two-leg ladders and the effect of the AFM inter¬ 
ladder exchange interaction on the low-temperature be¬ 
havior of these type of ladders. For an isolated spin-1/2 
FM two-leg ladder, a gapped behavior was observed at 
zero temperature. In principle, the ground state of an 
isolated spin-1/2 FM two-leg ladder is in the saturated 
FM phase. 

As soon as the inter-ladder AFM exchange interaction 
is added, a first-order commensurate-incommensurate 
quantum phase transition occurs in the ground state 
phase diagram. Numerical results of magnetization and 
the specific heat showed that the gapped FM phase is 
replaced by a nonmagnetic phase. In the presence of a 
magnetic field, the ground state of the coupled ladders 
remains in this nonmagnetic phase up to a critical field 
he- The value of the mentioned critical field depends on 
the value of the inter-lader AFM interaction. Although 
in the case of Cl( 0 -mi) 2 (-Cl) 2 )], the compound 

consists of isolated two-leg ladder, but such an enhance¬ 
ment of coupled interaction may occur upon chemical 
substitution. 

Existence of the non-magnetic phase in this model is 
very similar to what can be probably observed in the 
coupled spin-1/2 two-leg ladders with FM leg and AFM 
rung exchange interactions in the presence of a magnetic 
field2^“— . In fact, two-leg ladders with FM leg and rung 
interactions consisting of weak AFM inter-ladder cou¬ 
pling can also be considered as the two-leg ladders with 
FM leg and weak AFM rung exchange interactions which 
are FM coupled. The observation of non-magnetic phase 
in these new type of coupled ladders is stimulating for 
the future studies. 
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